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This paper is the companion article to [j4nn. Probab. 39 (2011) 
779-856]. We consider a discrete elliptic equation on the d-dimensional 
lattice Z'' with random coefficients A of the simplest type: They are 
Ph . identically distributed and independent from edge to edge. On scales 

(-H ' large w.r.t. the lattice spacing (i.e., unity), the solution operator is 

known to behave like the solution operator of a (continuous) ellip- 
tic equation with constant deterministic coefficients. This symmetric 
"homogenized" matrix Ahom ~ ahomid is characterized by ^ • AhomS, = 
({(, + V0) • A(£_ + S/(l))) for any direction ^ € R''-, where the random 
field (the "corrector" ) is the unique solution of —V* • A{^ + V(j}) = 
►^ ' in Z'' such that (f>{0) — 0, V(/> is stationary and {S7(f>) = 0, (•) denoting 

^Q , the ensemble average (or expectation). 

f*^ ' In order to approximate the homogenized coefficients Ahom, the 

^\ I corrector problem is usually solved in a box Ql = [—L, L)"^ of size 2L 

^^ . with periodic boundary conditions, and the space averaged energy 

-^ ' on Ql defines an approximation Al of Ahom. Although the statistics 

f^ , is modified (independence is replaced by periodic correlations) and 

r\l ■ the ensemble average is replaced by a space average, the approxima- 

tion Al converges almost surely to Ahom as L t oo. In this paper, we 
give estimates on both errors. To be more precise, we do not consider 
periodic boundary conditions on a box of size 2L, but replace the 
elliptic operator by T~^ — V* • ^V with (typically) T -^ L^ , as stan- 
$^ , dard in the homogenization literature. We then replace the ensemble 

^_' average by a space average on Ql, and estimate the overall error on 

the homogenized coefficients in terms of L and T. 
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2 A. GLORIA AND F. OTTO 

1. Introduction. 

1.1. Motivation. In this article, we continue the analysis we began in [6] 
on stochastic homogenization of discrete elliptic equations. More precisely, 
we consider real functions u of the sites x in a d-dimensional Cartesian 
lattice Z . Every edge e of the lattice is endowed with a "conductivity" 
a(e) > 0. This defines a discrete elliptic differential operator —V* • AV via 

-V* ■ {AVu)ix) := Yl aie)iuix)-u{y)), 

y&Z'',\x-y\=l 

where the sum is over the 2d sites y which are connected by an edge e = [x,y] 
to the site x. It is sometimes more convenient to think in terms of the 
associated Dirichlet form, that is, 

^Vi;-^Vu:= Y^ v{x){-V* ■ {AVu){x)) 
xeZ'i 

= X](^(^) ~ v{y))a{e){u{x) - u{y)), 

e 

where the last sum is over all edges e and {x,y) denotes the two sites con- 
nected by e, that is, e = [x,y] = [y,x] (with the convention that an edge is 
not oriented). We assume the conductivities a to be uniformly elliptic in the 
sense of 

a < a{e) < (3 for all edges e 

for some fixed constants < a < /3 < oo. 

We are interested in random coefficients. To fix ideas, we consider the 
simplest situation possible: 

{a(e)}e are independently and identically distributed (i.i.d.). 

Hence, the statistics are described by a distribution on the finite interval 
[a, /3]. We'd like to see this discrete elliptic operator with random coefficients 
as a good model problem for continuum elliptic operators with random co- 
efficients of correlation length unity. 

Classical results in stochastic homogenization of linear elliptic equations 
(see [8] and [13] for the continuous case, and [10] and [9] for the discrete case) 
state that there exist homogeneous and deterministic coefficients ^hom such 
that the solution operator of the continuum differential operator — V- AhomV 
describes the large scale behavior of the solution operator of the discrete 
differential operator —V* • AV. As a by product of this homogenization 
result, one obtains a characterization of the homogenized coefficients ^hom- 
It is shown that for every direction .^ E M'^, there exists a unique scalar 
field (j) such that V0 is stationary [stationarity means that the fields V<?i>(-) 
and V(/)(- -|- z) have the same statistics for all shifts z G Z'^] and (V(/>) = 0, 
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solving the equation 



(1.1) -V*-{A{^ + V(f))) = in 



yd 



and normalized by (/)(0) = 0. As in periodic homogenization, the function 
Z"' 9 a; I— )• ^ • a; + (pix) can be seen as the ^-harmonic function which macro- 
scopically behaves as the affine function Z 9 rr i— )■ ^ • x. With this "correc- 
tor" (f>, the homogenized coefficients ^hom (which in general form a sym- 
metric matrix and for our simple statistics in fact a multiple of the identity: 
^hom = flhomld) cau be characterized as follows: 

(1.2) c-Aome=((e+v</<)-A(e+v</))). 

Since the scalar field (^ + V(^) • A(^ + V(^) is stationary, it does not matter (in 
terms of the distribution) at which site x it is evaluated in the formula (1.2), 
so that we suppress the argument x in our notation. 

When one is interested in explicit values for Ai^orm one has to solve (1.1). 
Since this is not possible in practice, one has to make approximations. For 
a discussion of the literature on error estimates, in particular the pertinent 
work by Yurinskii [15] and Naddaf and Spencer [12], we refer to [6], Sec- 
tion 1.2. A standard approach used in practice consists in solving (1.1) in 
a box Ql = [—L, L)'^ n Z*^ with periodic boundary conditions 

(1.3) -V*-(A(e + V<AL,#))=0 in Ql, 
and replacing (1.2) by a space average 

(1.4) ^ • Al,#^ = Jq^ (e + V(t>L,#) ■ A{i + V(/.L,#) dx. 

Such an approach is consistent in the sense that 

lim Al,# = Ahom 

I/— >-oo 

almost surely, as proved, for instance, in [1] for the continuous case, and 
in [2] for the discrete case. Numerical experiments tend to show that the use 
of periodic boundary conditions gives better results than other choices such 
as homogeneous Dirichlet boundary conditions, see [14]. 

An important question for practical purposes is to quantify the depen- 
dence of the error (|Ahom — ^l,#P)^ in terms of L. Let us give another 
interpretation of (1.3): This equation on Ql is equivalent to (1.1) on Z'^ with 
a modified conductivity matrix Al, that is the periodization of A\q^ on Z*^. 
Doing this, we have replaced independent coefficients A by Qi-periodically 
correlated coefficients A. Since A and A are not jointly stationary (see Def- 
inition 4), it may be difficult to compare V(/) to V(J)l^^. To circumvent this 
difficulty, and following the route of [10, 13, 15] and [12], and as in [6], we 
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slightly depart from (1.3) by introducing a zero-order term in (1.1): 

(1.5) T~Vt-V*-(^(C + V(/)t)) = inZ"'. 

As for the periodization, this localizes the dependence of 4>t{z) upon A{z') 
to those points z' S U^ such that \z — z'\ < \/T (at first order). Yet, unlike 
the periodization, V(/>t and V(j) are jointly stationary. In terms of random 
walk interpretation, the lifetime of the random walker is of order T, and 
the distance to the origin of order yT. Hence, up to taking T ^ Lp' , in first 
approximation, the function 4>t\q^ only depends on the coefficients A{z) for 
z € Qli as it is the case for 4)l,#- 

We'd like to view 4't\Qt as a variant of </'l,# which is convenient for our 
analysis. We then define 

(1.6) e • Ar^Li = /" (e + V(/<t) • A{i + V<t>T)r]L dx, 

where tjl is a smooth mask with unit mass and support Ql. The aim of this 
paper is to determine the scaling of the error (|Ahom ~ At^lI'^)^'^ in terms 
of L and T. Eventually this will allow us to make a reasonable choice for T 
and L at fixed computational complexity. 

1.2. Informal statement of the results. When approximating vlhom by At^l, 
we make two types of errors: A "systematic error" and a "random error." 
In particular, as shown in [6], 

((e • AomC - e • AtM?) = (e • (Aom - At)0^ + var[AT,L]. 
The first term is the square of the systematic error (see [6], (1.10)) 

Errorsys(T) := m + V0t) • A{C + V0r)) - ((^ + V<A) • A{^ + V0))| 
(1.7) 

= {{V<pT - V(p) ■ A{V(j)T - V(/))). 

It measures the fact that the coefficient a(e) at bond e does (up to exponen- 
tially small terms) not influence (pTix) if |a; — e| ^ vT. This error vanishes 
for T = L^ t oo. The second term is the square of the random error, 

/ (e + V0r)-^(e + V</)T)r?Ldx . 

It measures the fluctuations of the energy density. This error vanishes as 
Ltoo. 

In [6], Theorem 1, we have proved that 



(1.8) Errorrand(r, L) = var 



(1.9) var 



(e + V(Ar)-A(e + V0T)r/LC?a; 



^/^^ id = 2, L^Mn^T, 



rs-/ 



d>2, L^'^/^, 



for some q depending only on a, /3, where "<" stands for "<" up to a mul- 
tiplicative constant depending only on a,/3 and d. We have also identified 
the systematic error in the limit of vanishing conductivity contrast, that is. 
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1 — j3/a <^ 1, and found 



Errorsys(r) ~ < 



(1.10) Errorsys(T) < 



d = 2, 


L-Hn'^L 


2<d<7, 


L-^/2, 


d = 8, 


L-^lnL, 


d>8, 


L-\ 



f d = 2, r-\ 

d = A, T-^lnT, 
d>4, T-2, 

where "~" means that both terms have the same scahng (in T). In this 
paper, we shall actually prove that for general a and (3 (see Theorem 1) 

d = 2, T-Mn^r, 

d = 3, r-3/2, 

d>4, T^^ 

where there is a logarithmic correction for d = 2 when compared to the 
vanishing conductivity asymptotics. 

Assuming that (J)t can be well approximated on domains of size L if we 
choose T ^ L^, the combination of (1-10) and (1.9) yields 



Ahom-AT,Lp)'/'< 



Hence, the numerical strategy converges at the rate of the central limit 
theorem for 2 < d < 8 (up to logarithmic corrections for d=2 and d = 8). 

Up to dimension 4, the systematic error for T '^ LP' scales as the square of 
the random error. In particular, this leaves room for the choice T. If we take 
T ^ L, then the systematic error is of the same order as the random error. 
What we have gained is that 4>t can now be well-approximated on domains 
of size R ~ \rr ^ vL, and not only L. Note also that the random error is 
unchanged if instead of taking the average of one realization of (f)T on Ql 
(with the mask /ix,) we take the empirical average of the averages of N 
independent realizations of (pT on a domain Q^/^i/d (with the according 
mask f^L/j^i/d). Hence, since <J)t can be well-approximated on domains of 
size R > vL, considering A^ = 1 realization of (pT approximated on Qi or 
N = vL'^ independent realizations of (pT approximated on Qrj^ yields the 
same scaling for the error between the homogenized coefficients and their 
approximations. Since the computational cost of solving a linear problem is 
superlinear in the number of unknowns, it seems best to choose N as large 
as possible, and therefore taking N = vL'^ seems a reasonable strategy at 
first order. Yet, we do not make precise in this paper the relation between R 
and vL in terms of absolute values (we only consider the scaling), which 
may make the optimal choice for A^ more subtle in practice than this general 
principle. A complete numerical analysis of the numerical method (including 
the influence of R and the optimization of A) will be presented in [4]. 



(1.11) ^.At^={^-AO- I 7TT^^der,iX). 
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We conclude this introduction by mentioning the very recent contribu- 
tion [11] by Mourrat. The equation under investigation is the same as above, 
namely a discrete elliptic equation on Z with i.i.d. coefficients. The object 
under study is the spectral measure associated with the generator of the 
environment viewed by the particle. Without entering into details, there ex- 
ists some nonnegative measure e^ associated with the elliptic operator and 
direction ,^ G M , such that the homogenized coefficient is given by 

e-Aome=(e-^o- / Ue.ix). 

Jr+ ^ 
As recalled in [11], we also have 

X-2/T 
{l/T + \f 
In particular, the systematic error can be written as 

i-Ari-i- Aon.e = ^ /^^ x{l/T + XY '^^»(^)' 

so that information on the scaling of the systematic error in terms of T 
yields information on the spectral behavior and conversely. The interplay 
between the strategy used in the present paper and the spectral measure is 
further investigated by Mourrat and the first author in [5]. In what follows, 
we do not make use of the spectral measure, which makes our approach 
self-contained. 

The article is organized as follows: In Section 2, we introduce the general 
framework and state the main results of this paper, that is, the systematic 
error actually scales as in (1.10). The last two sections are dedicated to its 
proof. 

Throughout the paper, we make use of the following notation: 

• d > 2 is the dimension; 

• J^d dx denotes the sum over x G Z'^, and jj^ dx denotes the sum over 
X G Z"' such that x £D, D subset of M'^; 

• (•) is the ensemble average, or equivalently the expectation in the under- 
lying probability space; 

• var[-] is the variance associated with the ensemble average; 

• cov[-; •] is the covariance associated with the ensemble average; 

• < and > stand for < and > up to a multiplicative constant which only 
depends on the dimension d and the constants a, /? (see Definition 1 below) 
if not otherwise stated; 

• when both < and > hold, we simply write ~; 

• we use ^ instead of > when the multiplicative constant is (much) larger 
than 1; 

• (ei, . . . , e^) denotes the canonical basis of Z . 
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2. Main result. 

2.1. General framework. 

Definition 1. We say that a is a conductivity function if there exist 
< a < /? < oo such that for every edge e of Z"^, one has a(e) G [a,/3]. We 
denote by Aa/s the set of such conductivity functions. 

Definition 2. The ehiptic operator L : Lf^^{Z'^) -!■ Lf^^{Z'^) , u^Lu 
associated with a conductivity function a G Aaji is defined for ah x G Z"^ by 

(2.1) {Lu){x) = -V* ■A{x)Vu{x), 

where 

u{x) — u{x — ei) 

V*u{x) :-- 



\/u{x 
and 



u{x + ei) — u{x) 
,u{x + e^) — u{x). 



.u{x) — u{x — Bd). 



A{x) :=diag[a(ei),...,a(ed)], 
ei = [x,x + ei],...,ed = [x,x + ed]. 

We now turn to the definition of the statistics of the conductivity function. 

Definition 3. A conductivity function is said to be independent and iden- 
tically distributed (i.i.d.) if the coefficients a(e) are i.i.d. random variables. 

Definition 4. The conductivity matrix A is obviously stationary in the 
sense that for all z G Z*^, A{- + z) and A{-) have the same statistics, so that 
for all X, z G Z"^, 

{A{x + z)) = {A[x)). 

Therefore, any translation invariant function of A, such as the modified cor- 
rector (pT (see Lemma 2), is jointly stationary with A. In particular, not 
only are (J)t and its gradient V(^t stationary, but also any function of A, <J)t 
and V(/t- a useful such example is the energy density {^ + Vcpx) • A{S, + V(/>t) , 
which is stationary by joint stationarity of A and ^(pT- 

Another translation invariant function of A is the Green functions Gt of 
Definition 6. In this case, stationarity means that Gr(- + 2;, ■ + z) has the same 
statistics as Gt{-, •) for all z G Z'^, so that in particular, for all x,y,z (i Z"^, 

{Gt{x + z,y + z)) = {GT{x,y)). 

Lemma 1 (Corrector; [10], Theorem 3). Let a G Aai3 be an i.i.d. con- 
ductivity function, then for all ^ G M"^, there exists a unique random func- 
tion (^ : Z'^ — 7- M which satisfies the corrector equation 

(2.2) -V*-A{x){^ + V(p{x)) = mZ'^, 
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and such that (p{0) = 0, Vcp is stationary and (Vcp) = 0. In addition, 

{\^<p\')<\e- 

We also define an "approximation" of the corrector as follows. 

Lemma 2 (Approximate corrector; [10], Proof of Theorem 3). Let a G Aa/s 
be an i.i.d. conductivity function, then for allT > and ^ G M'^, there exists 
a unique stationary random function (/ly : Z'^ — )■ M which satisfies the "ap- 
proximate" corrector equation 

(2.3) T-^^t{x)-V* ■A{x){C + V^t{x)) = inU^, 

and such that {(J)t) = 0. In addition, r^^((/'|.) + (|V0t|^) < |CP- 

Definition 5 (Homogenized coefficients). Let a S AafB be an i.i.d. con- 
ductivity function and let ^ € M and cf) be as in Lemma 1. We define the 
homogenized d x d-matrix ^hom as 

(2.4) e • ^home = ((^ + V0) • A(e + vct^m). 

Note that (2.4) fully characterizes ^hom since ^hom is a symmetric matrix 
(it is actually of the form Ohomid for an i.i.d. conductivity function). 

2.2. Statement of the main results. The main result of the article is the 
following estimate of the systematic error introduced in Section 1 . 

Theorem 1. Let a £ Aap he an i.i.d. conductivity function, and let <J)t 
denote the approximate corrector associated with the conductivity function a 
and direction ^ G W^ , \^\ = 1. We then define for all T ^ 1 the symmetric 
matrix At characterized by 

(2.5) ^^T^=((e + V</.T)•^(e + V</>T)>. 

Then, there exists an exponent q> depending only on a, f3 such that 
d = 2: |Aom-^T|<T-'(ln^)^ 

d = 3: |Aom-^T|<T^'/', 
(2.6) 

a! = 4.- |Ahom 

d>A: |Ahom-^r|<T' 



d = A: IAom-Arl<r-^lnr, 



As a by-product of the proof of Theorem 1 , we obtain the following corol- 
lary. 

Corollary 1. Let a G Aai3 be an i.i.d. conductivity function, d > 2, 
T > 0, and let cpx and <j) denote the approximate corrector and stationary 
corrector (see [6], Corollary 1) associated with the conductivity function a 
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and direction ^ G M , |^| = 1, respectively. Then 

(2.7) r-i((</.T-<A)') + (|V(/.T-V(A|2)<<^ d = 4, T-HnT, 

U>4, r-2. 

/n particular, 

lim (((0T - <^)') + (|V,/.T - V,^p)) = 0. 

This corollary gives a full characterization of the convergence of the reg- 
ularized corrector to the exact corrector for d> 2. 

Remark 1. Note that the definition (2.5) of At does not include the 
zero-order term T~^{(p'^), so that ^ • AtC does not coincide with the energy 
associated with the equation. Surprisingly, the addition of the zero-order 
term in the definition of At would make the estimate (2.6) saturate at T~^ 
for d>2. 

Remark 2. For d = 2, although we lose control of (pT we may still quan- 
tify the rate of convergence of V<^t to \^(p, the gradient of the corrector of 
Definition 1. In particular, (2.7) is replaced by 

{\V^T-V^\^)<T-hn'^T 

for some q> depending only on a,/3. 

2.3. Auxiliary lemmas. In order to prove Theorem 1 and Corollary 1, we 
need three auxiliary lemmas in addition to the results of [6] : The first one is 
a covariance estimate very similar to the variance estimate in [6], Lemma 2.3, 
the next one is a refined version of the decay estimates of [6], Lemma 2.8, 
whereas the last one is a generalization of the convolution estimate of [6], 
Lemma 2.10. 

Lemma 3 (Covariance estimate). Let a = {aj}jgp^ be a sequence of i.i.d. 
random variables with range [a,/3]. Let X and Y be two Borel measurable 
functions of a ^ R^ (i.e., measurable w.r.t. the smallest a-algebra on M'*^ for 
which all coordinate functions M^ 9 a i— )■ a^ E M are Borel measurable, cf. [7], 
Definition 14-4)- 

Then we have 



oo 



(2.8) cov[X;y]< V(sup 



i=l 



dX 



doi 



2\ 1/2 



sup 



dY 



doi 



2\ 1/2 



var[ai]. 



dZ\ 



where sup^. |^| denotes the supremum of the modulus of the ith partial 

derivative 

dZ 

— — [ai,. . . ,ai_i,aj,aj+i,. . .) 

oat 

of Z with respect to the variable ai G [a, (3], for Z = X,Y . 
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The proof of this lemma is standard. As for [6], Lemma 2.3, it rehes on 
a martingale difference decomposition. 

We define discrete Green's functions in the following definition. 

Definition 6 (Discrete Green's function). Let d>2. For ah T > 0, the 
Green function Gt -Aafs x Z'^ x Z'^ — ;• Z'^, {a,x,y) i— ;• GT{x,y;a) associated 
with the conductivity function a is defined for all y G Z and a G Aap as the 
unique solution G{-,y;a) G L^(Z'^) to 



T GT{x,y;a)v{x)dx + j Vv{x) ■ A{x)VxGT{x,y;a)dx = v{y) 
(2.9) 

where A is as in (2.1). 

Throughout this paper, when no confusion occurs, we use the short- 
hand notation Grix^y) for GTix,y;a). We need a decay of the Green func- 
tion Gxix^y) and its (discrete) gradient S/xGT{x,y) in \x — y\^ 1 that is 
uniform in a but nevertheless coincides (in terms of scaling) with the decay 
of the constant- coefficient Green function. The constant-coefficient Green 
function in the continuous case is known to decay as 

\x — y\ ~ expl — const. ^^ | for d > 2 and 

V Vt J 



— X — z/ \ 

In--^^ r 1 expl —const.-' — ^=^ for d = 2\ 

" y\) V Vt J 

its gradient decays as the first derivative of these expressions. Note the cross- 
over of the decay at distances |x — y| of the order of the intrinsic length scale 
vT ^ 1 from algebraic (or logarithmic in case of d = 2) to exponential. 

In the class of a-uniform estimates, these decay properties survive as point- 
wise in {x^y) estimates on the level of the discrete Green function Gxix^y) 
itself, but only as averaged estimates on the level of its discrete gradi- 
ent VxCrix^y). More precisely, VxGT{x,y) has to be averaged in x on 
dyadic annuli centered at x = y. It will be important that the average can 
be (at least slightly) stronger than a square average (see [6], Lemma 2.9). 
On the other hand, we do not need the exponential decay: Super algebraic 
decay is sufficient for our purposes. 

Lemma 4 (Pointwise decay estimate on Gt)- Let a G Aap, and Gt he the 
associated Green function. For d> 2, we have for all k> 0, and all x, y G Z"^ 

(2.10) GT(x,2/)<(l + |x-y|)2-'^min|l/'''"^'' 
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where the constant in '% " depends on k. For d = 2, we have for all k > 



(2.11) GTix,y)<{ 



In 



l + \x-y\ 



\x-y\ 



Vf 



for \x — y\ -^ vf 
for \x-y\>VT 



where the constant in "< " depends on k . 

Finally, for the proof of Theorem 1, we need to know that also the con- 
volution of the gradient of the Green's function with itself decays at the 
optimal rate, that is, with the following lemma. 

Lemma 5 (Convolution estimate). Let hT^gr'-'^ ^"^^ satisfy the fol- 
lowing properties. 

Assumptions on hnp [estimate of \^xGT{y + z,y)W: For all R^l and 

r>o, 



(2.12) 



d>2: 



R<\z\<2R 



hT{zfdz<R 



2-d 



(2.13) d = 2: 
and for i? ~ 1 

(2.14) d>2: 



I hT{zfdz<m:in{l,VTR-^f, 

JR<\z\<2R 



hT{zfdz<l. 
<\z\<R 

Assumptions on gx [estimate of Gxiy + z^y)]: For d > 2, and for all 

zeZ'^, 

-3- 



(2.15) 



gT{z) = {l + \z\y~''mm{l, 



and for d = 2, 

(2.16) gT{z) = { 

Then we have 



In I ^^ I for \z\ < VT 



1 —F= for \z\ > VT 

^ \Vf; 



(2.17) / griz) / hT{w)hT{z — w)dwdz 



(d = 2, T, 

d = 3, VT, 

d = i, InT, 

U>4, 1. 



3. Proof of the main results. Throughout this section, we let ^ G M"^ be 

such that 1^1 = 1. 
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3.1. Proof of Theorem 1. In view of (1.7), in order to estimate \At — 
^homl) we need to estimate how close the modified corrector (J)t is to the 
original corrector (f) [in terms of (|V(^t — V<^P)]- Therefore, it is natural 
to introduce tpx = T"^-^ (the prefactor T^ is such that 'ipT is properly 
renormalized in the limit T f oo at least for large d) . Considering ipj- is also 
convenient since for d = 2, the corrector (j) is not known to be stationary (only 
its gradient is known to be stationary) so that working with the modified 
correctors (pT, which are known to be stationary, avoids technical subtleties. 
In fact, we opt for a dyadically discrete version of ipT defined via 

(3.1) iPt-=T{(I)2T-4't)- 

This discrete version has the technical advantage that we do not have to 
think about the differentiability of <J)t in T. Moreover, its dyadic nature is 
in line with the dyadic decomposition of the T-axis according to 

oo 

(3.2) \At — ^homl < / ^ |^2'T ~ ^2'+itI 

forced upon us in the case of d = 2. In order to get (3.2), we used the fact 
that 

(3.3) lim AT = Ai,ora, 

T— !>oo 

which is proved in [6], Proof of Theorem 1, Step 8. We shall also use that ipT 
solves 

(3.4) T-^i^T - V* • AVijT = \>p2T- 

We split the proof in eight steps. 
Step 1. Derivation of 

(3.5) 1^ • {A2T - AT)i\ < r-2|(0rV'T)| + -1^\{<P2T^T)\. 

Although this could be directly inferred from the spectral formula (1-11) 
for At, we give an elementary argument relying only on the corrector equa- 
tion. We recall the following consequence of (2.3) which is proved in [6], 
Proof of Theorem 1, Step 8: 

(3.6) r-i(0TX) + ((e + V(AT)-AVx)=O 

for every field x : Z"' — t- M that is jointly stationary with A and such that 
(x^) < oo. From formally differentiating the definition (2.5) of At w.r.t. T 
and using (3.6) for x = ~§f'^ ^^ obtain 

We claim that the corresponding discrete-in-T version reads 

(3.7) e • {A2T - At)^ = -T-\{i,TM + U^T<P2T))- 



OPTIMAL ERROR IN STOCHASTIC HOMOGENIZATION 



13 



Indeed, by definition of At, by expanding tlie square, by symmetry of A, 
by definition of ipT, and (3.6), we have 

e • {A2T - At)^ 

= ((e + V(A2t) • A{^ + Vct>2T)) - ((e + V^t) ■ A{^ + Vct>T)) 

= ((V(^2T - V(/.t) • A{^ + V<t)2T)) + ((V(/>2T - V0t) • ^(C + V0t)) 

('i^ t-\{v^Pt • ^(e + V02t)) + (VV'T • ^(e + V(/.t))) 
^=^ -r-i((2r)-i(^T</'2T) +t-1(^t</'t)). 

In the next four steps, we focus on the first term of the r.h.s. of (3.5). The 
second term will be dealt with the same way in Step 7. 
Step 2. Proof of 



(3.8) K</'TV'r)|<j;(sup 

g \a(e) 



dMo) 



dale) 



2\ 1/2 / 



(sup 

\a(e) 



dMO) 



da{e) 



2\ 1/2 



where the sum runs over the edges e, and proof of the representation formulas 
9<At(0) 



(3.9) 



(3.10) 



da{e) 

gV'T(O) 
5a(e) 



-(ei + V,</.T(^))V,,GT(^,0), 

-\ I GT{QM{ii+'^t<t>2T{z))yzfi2T{z,w)dw, 

^ .h.d 



where the edge is e = [z, z + e^] . 

Due to [6], Lemma 2.6, the functions (pT and ipT are measurable with 
respect to the coefficients a. Hence, (3.8) is a consequence of the covariance 
estimate of Lemma 3: Since {(px) = {^Pt) = 0, 

{(pTlpr) = {{4>T - {(t)T)){i^T - {iPt))) 

= cov[(J)t;'>Pt]- 

Formula (3.9) is identical to [6], Lemma 2.4, (2.12). To prove (3.10), we 
first make use of the Green representation formula for the solution to (3.4): 



(3.11) 



iPt{x) 



Gt{x, w)4)2t{w) dw 



for all X G Z*^. Since a(e) i— )• (j)T{-',a{e)) and a(e) i— )• (/)2T{']o,{e)) are contin- 
uously differentiable by [6], Lemma 2.4, we deduce by formula (3.1) that 
a(e) I— 7- ipTi'^cLi^)) is also continuously differentiable. Using then the formu- 
las [6], Lemma 2.5, (2.15), and [6], Lemma 2.4, (2.12), for the derivatives 
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of Gt and (I)t with respect to a(e), and the fact that Gt G L^(Z ) (see [6], 
Corollary 2.2), we may switch the order of the differentiation and the inte- 
gration to obtain for all x G Z'^ 

dipxix) 1 f dGxixjw) 



4>2t{w) dw 

1 r 
+ - / Gt{x,w) 1^''/r^ dw 

2 Jzd 



da{e) 2 J^d 9a(e) 

d(j)2T{w) 



da{e) 

f-i-tn^ [6],(2.12) and (2.15) If \V7 n ( \A. ( \ ^ 

(3.12) = -- / V zfiT{x,z)V zfiT[z,w)cl)2Tiw) dw 

-n GT{x,w){ii+Vi(j)2T{z))Vzfi2T{z,w)dw 

2 Jl'i 

^='^-V,,Gt(x,z)V,^t(^) 

- T, I GTix,w){^i + \7i(l)2T{z))^z,G2T{z,w)dw, 
2 Jzd 

which is (3.10) taking x = 0. 

From now on in the proof, we let qt be defined as in Lemma 5 (i.e., gx 
decays as the Green function Gt)- 

Step 3. In this step, we shall prove that 

(3.13) |(0tV't)|</: + AA, 
where 

C := / ((1 + |V</.r(z)P)|V,GT(z, 0)\Y' 

(3.14) ^' 

^ ^ 2\ 1/2 



(3.15) 



gT{w){l + \V<P2T{z)\)\V,G2T{z,w)\dw] ) dz, 

and A/' = A/'i +A/'2, 

M:= / ((1 + |V<At(^)P)|V,Gt(^,0)|2)1/2 

x(|V7AT(^)P|V,GT(z,0)|2)l/2d^, 

Ar2:=Mr) / {il + \VMz)\')\VzGTiz,0)\Y' 

x{{l + \V^2T{z)f)N,GT{z,0)\'f^dz 



(3.16) 
with 
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The term £ is a linear error: It is of the same type as for the analysis 
in the limit of vanishing ellipticity contrast (see [6], the Appendix). On 
the contrary, the term N is nonlinear and does not appear in the limit of 
vanishing ellipticity contrast. As we shall prove, it is of lower order. The 
terms C and A/i in estimate (3.13) would be direct consequences of (3.8), 
and (3.9) and (3.10), disregarding the suprema in a(e) in (3.8). Taking the 
suprema in a(e) into account actually brings the second nonlinear term A/'2, 
which turns out to be of lower order than A/i . 

According to [6], Lemma 2.4, (2.13), we have for (3.9) 



(3.17) 



sup 

a{e) 



OT 



(0) 



da{e) 



{1 + \V^4>t{z)\)\^zGt{z, 



It remains to deal with (3.10). Using the pointwise decay of Gt in Lemma 4 
combined with the susceptibility estimates [6], Lemma 2.4, (2.14), and [6], 
Lemma 2.5, (2.16), of V(^r and VGt w.r.t. a(e), we obtain 



(3.18) 



sup 

a(e) 



GT{^,w){i, + V42T{z))Vzfi2T{z,w)dw 



< / gTHi'^ + \V^(|)2Tiz)\)\V,G2Tiz,w)\dw, 



which together with (3.17) gives the linear term C 

To treat the first term of the r.h.s. of (3.10), we need to deal with the 
supremum of |Vj^r(z)| over a(e). We appeal to (3.12) that we rewrite in 
the form 



dipT{x) 
da{e) 



-Vi^PTiz)GTix,e) 



7:(Ci + Vi(/>2T(^)) / GT{x,w)GTie,w)dw, 

2 Jzd 



where GT{x,e) := Gt{x,z + e^) — Gxixjz) and GTie,w) := Gt{z + ei,w) 
Gt{z^w). Hence, 



(3.19) 



da{e) 



-V,M^)GT{e,e) 



7;i^i + ^i4>2T{z)) / GT{e,w)G2T{e,w)dw, 
2 Jzd 



where GT{e,e) := GT{z + ei,z + ei) + GT{z,z)-GT{z + ei,z)-GT{z,z + ei). 
On the one hand, the uniform bound [6], Corollary 2.3, on VGt yields 
|Gr(e, e)| < 1. On the other hand, as we shall argue, the integrability of VG^ 
and VG2T from [6], Lemma 2.9 (combined with the uniform bound [6], 
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Corollary 2.3, on gradients) implies 

(3.20) ^^GT(e,^)G2T(e,^)ci«;<^,(r) = |^ = ^, InT, 

Hence, if we regard (3.19) as an ordinary differential equation for VjV'r(-2) 
in the variable a(e), we obtain 

(3.21) SUp\ViMz)\<\'^^Mz)\+^^d{T){l + \V^(j)2Tiz)\) 
a{e) 

since a(e) lies in a bounded domain [a,/3], and sup^j/g) |Vj(/)2t(-z)| ^ 1 + 
\^i(p2T{z)\ according to [6], Lemma 2.4, (2.14), with 2T instead of T. Note 
that (3.17), (3.21) and sup„(e) \Vz,Gt{z,0)\ < \VzfiT{z,0)\ give the nonlin- 
ear terms A/i and A/2 • 

We now give the argument for (3.20). We first use the Cauchy-Schwarz 
inequality 

Grie, w;)G2T(e, w) dw 
<( [ GTie,wfdw] ( [ G2T{e,wfdw 



1/2 / /. \ 1/2 

<(/ \V,GT{z,w)\^dw] / \V,G2T(.Z,w)\^dw 



and then make a decomposition of Z into the ball of radius iZ ~ 1, and 
dyadic annuli {w : 2^R <\z — w\ < 2*+^i?} for i E N. On the ball of radius R, 
we use the uniform estimate of [6], Corollary 2.3, on VGt, whereas on the 
dyadic annuli we appeal to the decay estimate in [6], Lemma 2.9, for the 
gradient of the Green function, which requires R to be sufficiently large 
although still of order 1. Both terms in the r.h.s. scale the same way and we 
only treat the first one: 

\V;,GT{z,w)\'^dw 
= / \V,GT{z,w)?dw + y^ \V,GT{z,w)?dw 

J\z-w\<R j^Q J2»_R<|2-io|<2»+i_R 

oo 

< 1 + ^(2*)'='+2(i-'^) min{l, VT{2'R)-^f 



</"d(r), 



using [6], Corollary 2.3, and [6], Lemma 2.9, for k = 2, respectively. This 
concludes Step 3. 
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Step 4. Suboptimal estimate of the nonlinear term A/": 

^ (d = 2, VTln'^T, 

(3.22) Ml < (VVt • AV^Jjt)^/^ {d = 3, InT, 

[d>3, 1, 

(3.23) AA2</Xrf(T)^ 

where q is a generic exponent which only depends on a,/3. We first deal 
with A/i, and begin with the second factor of the r.h.s. of (3.15). The point- 
wise estimate (2.10) of Lemma 4 for d> 2 on the Green function gives the 
suboptimal pointwise estimate on the gradient of the Green function 

d 

(3.24) \VGt{z, 0)1 < Gt{z, 0) + J^ Gt{z + e„ 0) < (1 + \z\f-''. 

This estimate coincides for d = 2 with the uniform bound of [6] , Corol- 
lary 2.3. The coercivity of A thus yields 

{\VGT{z,0)\^\VMz)\y^' 

<il + \z\f~''{VMz) ■ A{z)VMz))'^' 
= (1 + |z|)2-'^(VVt-^V^t)'/' 
by joint stationarity of VipT and A. Hence, (3.15) turns into 

AAi<(VVt-^VV^t)'/' / {l + \z\f-''{{l + \VM^)\^)NGT{z,0)\y^^dz. 

We then let p > 2 be a Meyers' exponent as in [6], Lemma 2.9 and use 
Holder's inequality in probability with exponents {p/{p — 2),p/2), the sta- 
tionarity of V(j)T, the fact that the gradient of (pT is estimated by (pT as in 
(3.24), and the bounds on the stochastic moments oi (pT in [6], Proposition 1, 

X / (l + |z|)2-'^(l + |V</.r(z)|2P/(P-2)^(p-2)/(2p)^|^^^^^^Q^|p^l/p^^ 

(3.25) X / (l + |z|)2-'^(|VGT(z,0)r)i/Pdz 

< {ViPt ■ ^VVt)'/'(1 + |,^^| V(p-2)^(p-2)/(2p) 

X / [l + lzlf-^ilVGrizMn'^^dz 



<UTY{yi'T-Ay^T)^'^ / (l + |z|)2-'^(|VGT(z,0)nVP(iz, 



for some generic q depending only on a,/3. Holder's inequality with expo- 
nents {p,p/{p— 1)) in Z,'^, combined with the same dyadic decomposition 
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of Z as for the proof of (3.20) (and the uniform bound on VGt from [6], 
Corollary 2.3) yields 

1/p 

VGT{zM''dz)\ 

i=0 "^ J2iR<\z\<2i+^R 



^1 + Ef(/ WGTizM'dz)) 

x( [ (l + |z|)(2-^)P/(P-^)ci^ 

\J2iR<\z\<2i+^R 



{p-l)/p 
l + |^|)(2-rfW(p-l)d^) 

/2»_R<|2|<2»+iH 

Using the optimal decay of VGt on dyadic annuli in L^ norm from [6], 
Lemma 2.9, with k = 2p, this turns into 

{l + \z\f-''{\VGT{z,0)n'^''dz 

! 



<l + 5;(^(2^i2)'^(2'i?)(^"'^)^min|l,|^| J 

x((2^i?)'^(2^i?)(2-«')p/(p-i))(P-i)/P 
= l + f;(2^i?)3-min{l,^f. 

i=0 "" ^ 

Recalling that i? ~ 1 , this implies 

/■ , f(i = 2, ^/T, 

/ il + \z\f-''{\VGT{z,0W)'^'dz<{d = 3, InT, 
•'^' [d>3, 1. 

Combined with (3.25) it proves (3.22). 

We now turn to N2- Proceeding as above to deal with the terms ^(pT 
and V(^2T in J^2 ■, we obtain as desired 



AA2<M,(T)/.,(T)25 / (|V,GT(^,0)r)2/Pdz 



yd 



using the same dyadic decomposition of Z, as for the proof of (3.20) to- 
gether with the higher integrability of gradients of [6], Lemma 2.9 and [6], 
Corollary 2.3. 

Step 5. Estimate of the linear term C: 

(d = 2, Tln^r, 
d = 3, VT, 
d = 4, InT, 
U>4, L 



(3.26) £ < <^ 
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We first treat the second factor of (3.14). We proceed as in Step 4 to deal with 
the expectation of the corrector term, and let p > 2 be a Meyers' exponent 
as in [6], Lemma 2.9. We obtain by Holder's inequality in probability with 
exponents {p/{p — 2),p,p) and the bounds on the stochastic moments of (J)t 
from [6], Proposition 1: 

grHil + \V(l)2Tiz)\)\V z,G2Tiz,w)\dw' 



9T{w)gT{w') 

I 

x{{l + \V(l,2T{z)\f\yzfi2T{z,w)\\V,fi2T{z,w)\)dwdw' 

gT{w)gT{w'){\Vzfi2T{z,w)\P)^'P 
x{\Vzfi2T{z,wT?'^''dwdw' 
<f^d{TrU9T{w){\VzG2T{z,wr)'/Pdw) . 
We thus have 

^<w(rr / {ii + \vMz)\')\VGTiz,o)\Y' 

X / gT{w){\VzG2T{z,w)n'/Pdwdz. 



Appealing once more to Holder's inequality in probability with exponents 
{p/{p — 2),p/2) and to [6], Proposition 1, this turns into 

^<f^d{Tf' [ 9t{w) [ {\VzG2T{z,w)n'/P{\VzGT{z,0W)'^''dzdw 



= lid{Tf / griw) / h2T{z-w)hT{z)dzdw, 

where, by stationarity, we have set 

hT{w) = l,\V^GT{wM''?'\ 

By the optimal decay estimate of VGt on dyadic annuli from [6], Lemma 2.9 
(and by the uniform bounds on VGt from [6], Corollary 2.3), and by defi- 
nition of gT^ we are in position to apply Lemma 5. Estimate (3.26) is thus 
proved. 
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Step 6. Proof of 

(3.27) (VVT-^VVT)<|(</'rV'T)|. 

Using (3.1), we rewrite (3.4) as 

(2r)" Vr - V* • AVVt = U2T - (2T)^ Vr 
(3.28) 

We now multiply (3.28) by ^l^T'■ 

(2r)~ Vr - (V* • AVtljT)^T = ^T^T- 

By integration by parts and joint stationarity of V'T; V^t and A (see [6], 
Proof of Theorem 1, Step 8, for details), this turns into 



We then conclude by the nonnegativity of the first term. 


Step 7. Proof of 




■ d = 2, 


Tln'^T, 


(3.29) \{<PtM\<< ji^; 


Vt, 

InT, 


(d>A, 


1, 


and 




(d = 2, 


Tln^T, 


(3.30) K</>2tV't)|<< 


d = 3, 

d=4:, 


VT, 

InT, 




.(i>4, 


1. 


From Steps 3, 4 and 5, and Young's inequality, we deduce that 




fd = 2, rin«r, 


K</'rV'r>|--(V^T-^VV'T)< < 


d = 3, Vf, 
d = A, InT, 






.d>4, 1. 



Combined with Step 6, this shows (3.29). 

For (3.30), we proceed exactly as for (3.29) in Steps 2-6. In particular, 
with obvious notation, we have 

\{<t^2T^T)\<N' + C', 

where 

, (d = 2, TlniT, 

M'-- ( VV^T • ^VV^t) <{d = 3, ln2 T, 

^ U>3, 1, 
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a < < 



and 

d = 2, TlnT, 

d = 3, VT, 

d = 4, InT, 

U>4, 1. 

We then conclude as above. 
Step 8. Proof of (2.6). 
Steps 1 and 7 yield 

\C-{AT-A2T)C\<T~^\{(I^T^I^T)\ + {2T^r^\{cl)2T4^T)\ 

(3.31) 





'd = 2 
d = 3 
d = 4 




TWT, 

Vt, 

InT, 




.d>4 




1. 


We finally appeal to the dyadic decomposition of the T-axis (3.2), which, 
combined with (3.31), turns into 


' d = 2, 

oo 
|^(AT-^hom)e|<>^]^-4' 

U>4, 




(2*r)-Mn'?(2^r), 
(2^r)-3/2, 
(2'r)-^in(2*r), 
(2^r)-2^ 


< . 


■d = 2, 
d = 3, 
d = 4. 


T- 
T- 


-iln^T, 

-3/2 

-2 InT, 




.d>4. 


T- 


-2 


This concludes the proof of the 


theo 


rem. 







3.2. Proof of Corollary 1. By Steps 6 and 7 in the proof of Theorem 1 
and by the definition (3.1) of ipx, we learn that 



(|V<^2T-V0rP 



(3.1) 



t-Wv^t\ 



(3.27) and (3.29) f d = 3, 
< I d = 4, 



2^-3/2 



d>A, 

In particular, \/(pT is a Cauchy sequence in L^ in probability. Hence, \/(J)t 
converges in L^ to its weak limit V(j), and by a dyadic decomposition of the 
T-axis the above estimate yields 

(d = 3, T-^n^ 
(|V0T-v,/.|2)< J d = 4| r-2inr, 
[d>4, r-2, 

which gives the second term of the l.h.s. of (2.7). 
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(3.29) and (3.30) f d = 3, 
< <'d = 4, 



Likewise, from Step 7 in the proof of Theorem 1, we learn that 

{{(j)2T - (prf) = T~^ {{(t)2T - (t>T)^T) 

<T-\{\ct)2T^T\) + {\(t^T^T\)) 

T-MnT, 
d>4, T-i, 

so that (l)T is a Cauchy sequence in L?' in probabihty and c/r converges in L^ 
to its weak hmit (p provided by [6], Corohary 1. In particular, by a dyadic 
decomposition of the T-axis the above estimate yields 

d = i T~HnT, 
d>4, T-\ 

which is the first term of the l.h.s. of (2.7). This concludes the proof of the 
corollary. 

4. Proof of the auxiliary lemmas. 

4.1. Proof of Lemma 3. Without loss of generality we may assume 



(4.1) 



oo 

^\sup 

a 



dX 



doi 



dY 



oo / 
>,^(sup 7^ ) < OO. 

Let Zn denote the expected value of Z conditioned on ai, . . . ,a„, that is 

Z„(ai,...,a„) := (Z|ai, . . . ,a.„). 
From [6], (5.2) and (5.3), in the proof of [6], Lemma 2.3, we learn that 

lim((Z-Z„)2)=0 

ntoo 

for Z = X,Zn = Xn and Z = y, Z„ = y^, respectively, so that, by the Cauchy- 
Schwarz inequality in probability, 

lim(X„) = (X), 

nfoo 

iim(y„) = (y), 

nfoo 

iim(x„y„) = (xy). 

nfoo 



Hence, 



(4.2) 



hm cov[X„;y„] = lim((X„y„) - (x„)(y„)) 

nfoo nfoo 



{XY)-{X){Y) 
cov[X;y]. 



OPTIMAL ERROR IN STOCHASTIC HOMOGENIZATION 23 

Note also that 

n 

(4.3) cov[X„; y„] = ^((X,y,) - (X.^iYi^i)) 

with the notation Xq = {X) and Yq = {Y), so that {Xn) = Xq and (Yn) = Yq. 
Inequality (2.8) then follows from (4.1), (4.2), (4.3), and 

2\ 1/2 



dX 



dtti 



2\ 1/2 / 

(sup 



dY 



dtti 



(4.4) (x,y,)-(x,_iy,_i)<(sup 

that we prove now. By our assumption that {ajjjgpj are i.i.d., we have 
Zj_i(ai,...,ai_i) = / Zi(ai,...,ai_i,a-')/?(da-'), 
{Xi{ai, . . . ,ai)Yi{ai, . . . ,ai)) 



Xi{ai,..., aj_i, a-)yj(ai, . . . ,aj_i,a-)/3((ia-) V 
where /? denotes the distribution of oi. Hence, we obtain 

{XiY,) - (Xi.iy^.i) 

Xj(ai, . . . , ai-i,a'i)Yi{ai, . . . ,ai_i,a-)/3((ia-) ) 

Xi(ai,...,aj_i,a-)/3(da-) / yj(ai, . . . ,aj_i,a'/)/3((ia-') \ 



:{Xi{ai, ..., aj_i,aj) - Xi(ai, . . . , aj_i, a-')) 



X (yj(ai, . . .,ai-i,a'i) - Yi{ai, . . . , ai-i,a'l))l5{da[)j3{da'l) 



< 



1 \ 

-(Xi(ai, . . . , Oj-i, a-) - Xi(ai, . . . , Oj-i, a-'))^/3(c^a-)/3(da'/) j 

-iYi{ai, . . . ,aj_i,a-) - Yi{ai,. . . ,ai-i,al)f I5{da[)l3{dal)\ 



1/2 



We then conclude the proof of (4.4) as in the proof of [6], Lemma 2.3. 

4.2. Proof of Lemma 4- We divide the proof in two main parts and deal 
with \z\ < vT and \z\ > yT separately. The proof relies on the Harnack 
inequality on graphs. We refer to Zhou [16] for Z , and to Delmotte [3] for 
other graphs. We recall here the easy part of Harnack's inequality (see [3], 
Proposition 5.3, or [16], Proof of Theorem 3.3, (3.11)). 
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Lemma 6 (Harnack's inequality). Let a £ Aa/s and i? » 1. //^r : Z*^ — ^ M+ 
satisfies 

(4.5) - V* • AVg{x) < 

in the annulus {R/2 < \x\ < AR} (i.e., g is a nonnegative suhsolution), then 

/ r \ 1/2 

(4.6) sup ff(x)<LR"M g{xfdx 

R<\x\<2R \ JR/2<\x\<AR 

Step 1. Proof of (2.10) for |x - y| < ^/T. 
Since Gt satisfies 

(4.7) - V^ • AV,,Gt{x, y) = -T-^Gt{x, y) < 
for |a; — y| ^1, one may apply Lemma 6. For R^ 1, we then have 

sup GT{x,y)<{R~'' Grix.yfdx 

x:R<\x-y\<2R \ J R/2<\x-y\<4R 

Combined with [6], Lemma 2.8, (2.21), for q = 2 (which is uniform in T > 
and y G Z ), this yields 

sup GTix,y)<R^-'', 

R<\x-y\<2R 

from which we deduce (2.10) for vT > |x — y| ^ 1. For |x — y| ~ 1, we 
appeal to [6], Proof of Lemma 2.8, (4.4), with i? ~ 1 and q = l, which yields 
sup\ ^_y\^f(GT{x,y) ^ 1 by the discrete L^ — L°° estimate. 

Step 2. Proof of (2.11) for |x - y| < y/T. 

Let A^ be a positive integer such that 2 ~ vT and 2~ vT S> 1. For all 
i G {1, . . . , A^}, we first show that 

({2-'^/Tr'' [ GT{x,yfdx^ 

\ J2-^-'^VT<\x~y\<2'i+^VT 

(4.8) 

< i ~ In 



i + 2-»Vr 

Estimate (4.8) follows from the triangle inequality and the BMO estimate 
of [6], Lemma 2.8 (2.20), provided we show that 

(4-9) GT{\^_y\^2-^+2^y<i, 

where Gth^^ \^2~i+^y/T\ denotes the average of Gxix^y) on the set {|x — 
y\ < 2^*''~^vT}- By the triangle inequality and the BMO estimate of [6], 
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Lemma 2.8, (2.20), we have 

^T^ {\x-y\<2-^+'^VT} 

- ^r{|a._y|<2-i+3^} 

+ 2' 



x-y\ <2-*+3VT}| 

1/2 



where C is a universal constant independent of i. Combined with the esti- 
mate for i = l 

^'^{|a:-?/l<4VT} ~ ^' 

which is a consequence of [6], Lemma 2.8, (2.22), this imphes (4.9) by in- 
duction. 

We are now in position to prove (2.11) for \x — y\ < vT. Since x i— )■ 
GT{x,y) satisfies 

-V* • AVxGTix,y) = -T-^GT{x,y) < 

in the annulus {x, 2~*~^\/T < \x — y\ < 2~*^^\/r}, Lemma 6 imphes 

sup GT{x,y) 

x:2-iVT<\x-y\<2-i+'>-VT 

<((2-'Vt)-2^ ^ _ _GT{x,yfdx^ 



<ln 



VT<\x-y\<2-'+'2VT 

T 



,1 + 2-*\/T, 

using (4.8) for 2~*vT » 1. For \x — y\ < R^ 1, we appeal to (4.9) and to 
the discrete L^ — L°° estimate 

GT{x,y) <i?^GT{|:c-j;|<R} <lnT. 

This completes the proof of (2.11) for \x — y\ < \/T. 
Step 3. Proof of (2.10) and (2.11) for \x - y\ > VT. 
Let R > \fT . Since Gt satisfies 

-V:. • AVxGrix.y) = -T-^Gt{x, y) < for \x-y\> 1, 

Lemma 6 implies 

sup GT{x,y)<(R~'' [ GT{x,yfdx^ 

x:R<\x-y\<2R \ J R/2<\x-y\<4R j 
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Combined with [6], Lemma 2.8, (2.23), for q = 2 and r = k, that is, 



f GT{x,yfdx<R''^^''-''^WTR' 

J Rl2<\x-v\<iR 



-l\k 
I R/2<\x-y\<4:R 

this yields the desired pointwise bound. 



4.3. Proof of Lemma 5. First note that by symmetry 

/ hT{z)hT{z — x)dz = / hT{z)hT{z — x) dz 

J\z\<\z-x\ J\z\>\z-x\ 

— n hT{z)hT{z — x)dz. 
Hence, it is enough to consider 



grix) / hT{z)hT{z — x)dzdx. 

' J\z\<\z-A 

In this proof, we essentially combine the pointwise decay of qt with the 
results of [6], Lemma 2.10, that we recall here for the reader's convenience 
(see [6], Proof of Lemma 2.10, Steps 1, 2 and 4): There exists i? ~ 1 such 
that for all R > R/2, 



hT{z)hT{z — x) dzdx 

jR<\x\<2RJ\z\<\z-x\ 
(4.10) 

<fd = 2, /22^ax{l,ln(Vri?-i)}, 
~ld>2, R^, 

(4.11) / _ / hT{z)hT{z-x)dzdx< ^'^ = ^' ^^^' 



\x\<4RJ\z\<\z-x\ L'^>2, 1. 

In view of (4.11) and (4.10), it will be convenient to make a dyadic de- 
composition of space. In order to also benefit from the decay of grix) for 
\x\'^ vT, we make the following decomposition of Z : 

(4.12) Z'^ = {\x\<2-^Vt} 

(4.13) U y {2'y/T<\x\<2'+^VT} 

(4.14) u|J{2*\/T<|x| <2^+1\/T}, 

where / is characterized by 2R < 2~^ \/T < 4R. 

For the integral over the r.h.s. of (4.12), we appeal to (4.11) and to the 
definitions (2.15) and (2.16) oi gT{x) for \x\ < 1: 

d = 2, In^T, 



(4.15) / grix) hT{z)hT{z-x)dzdx< 

J\x\<2-'Vt J\z\<\z-x\ la> z, i. 
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For the integral over (4.14), we use this time (4.10) for R > yT and the 
definitions (2.15) and (2.16) of gri^) for \x\ > vT, so that for all i G N we 
have 



gri^) / hT{z)hT{z — x)dzdx 

2»VT<|x|<2'+iVT J\z\<\z-x\ 

Summing this inequality on i € N then yields the estimate 



(4.16) 



gri-c) / hj'{z)hT{z — x)dzdx < VT 

VT<\x\ J\z\<\z-x\ 



A-d 



We now deal with the integral over the last part (4.13) of Z"^. To this 
aim, we combine (4.10) for R < vT with the definitions (2.15) and (2.16) 
of grix) for |x| < vT- In particular, for all i G {— /, . . . , —1}, we have 



grix) / hT{z)hT{z — x)dzdx 

2iVT<\x\<2i+'>-VT -'|^|<|2-2:| 



< 



d = 2, ln{2-'){2'VTf ln(2-*) ~ ^2(2*^)2^ 
d>2, {TVTf{2'VTf-'^ = {2'VT)^-'^. 

Summing this inequality over i G {— /, ■ • . , —1} and using that 2 ~ yT then 
yield 



I 



2-^VT<\x\<VT 

d = 2 



gxix) I hx{z)hT{z — x)dzdx 

\z\<\z-x\ 



(4.17) 



< 



1 + r ^ i^A' 



i=-i 



4-d 



d>2, l + \/r ^(2 



4-d\i 



2'r^VT 
< 



d = 2, 

d = 3, 

d = 4, 

l(i>4, 



VT, 

InT, 
1. 



The combination of (4.15), (4.16) and (4.17) finally proves (2.17). 
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